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Abstract. A microscopic model of a Josephson junction between two superconduct- 
ing plates is proposed and analysed. For this model, the nonequilibrium steady state of 
the total system is explicitly constructed and its properties are analysed. In particular, 
the Josephson current is rigorously computed as a function of the phase difference of 
the two plates and the typical properties of the Josephson current are recovered. 
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1. Introduction 

In 1962, Josephson |1| predicted a novel phenomenon in superconductivity, namely when 
two different superconductors were brought into close contact. Based on elementary 
quantum mechanics, he predicted the existence of a supercurrent with a peculiar current- 
voltage dependence. Namely, he argued that there would emerge a current of Cooper 
pairs which is proportional to the sine of the phase difference of the order parameters of 
both superconductors. The success of this prediction was immediate when indeed this 
phenomenon was experimentally observed already one year later 0. It counts as one of 
the greatest successes of quantum mechanics in physics and you will find a chapter on the 
Josephson effects in almost every textbook on superconductivity. During the following 
decades, the increase of knowledge on this subject in theoretical solid state physics 
has been tremendous and applications of Josephson junctions in electronic devices have 
been developed pi. Progress in conceiving a microscopic theory for the Josephson 
effect in rigorous quantum statistical mechanics was made when Sewell obtained the 
Josephson and Meissner effects in an model independent approach from the assumption 
of offdiagonal long range order and local gauge covariance 

Here, we are concerned with a substantive and rigorous derivation of the Josephson 
current in a concrete microscopic quantum statistical model. The model consists of two 
twodimensional superconducting plates which have a common onedimensional contact 
surface through which Cooper pairs can tunnel, and as such it induces a current. We use 
recently developed ideas of nonequilibrium statistical mechanics to conceive a reasonable 
framework for this model in which the Josephson current can be computed explicitly. 
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The model at hand is based on the mean field version, sometimes called the 
Anderson version, of the strong coupling BCS model. This model may lack many 
physical features, it is nevertheless an interesting model as it allows rigorous and explicit 
computations leading to nontrivial results. 

Due to the phase difference between the superconductors, the total system is not in 
an equilibrium state and we have to find a reasonable description of the nonequilibrium 
steady state of the system. In spite of the great success of statistical mechanics for 
systems in equilibrium, still not very much is understood in the nonequilibrium situation, 
especially for quantum systems. One of the most important points in this note is thus 
the construction of a good nonequilibrium steady state for this model (Section I2.2|l . 
This nonequilibrium steady state is constructed by selecting the state which fulfils a 
list of natural conditions in this setup. This construction is then compared with other 
treatments found in the hterature |S1IE]- 

After developing our framework for this model, we calculate the currents in the 
system (Sectional). First, we find that there is a Cooper pair current with the typical 
phase dependence of the traditional Josephson currents (see Figure Ej). Secondly, this 
model can be compared with other nonequilibrium quantum statistical models for 
entropy production [HllZllHlini and we also compute the heat fiuxes between the two 
superconductors. In spite of the fact that there are nontrivial particle currents in this 
system, we get that there are no heat fiuxes or entropy production here. 



2. The microscopic model 



The model we propose consists of two superconducting plates, I and //. At one side, 
the distance between these superconductors is made very small, such that Cooper pairs 
can tunnel through this barrier, see Figure ^ i.e., the superconducting plates have a 
common contact surface. 
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Figure 1. Two superconductors with a contact surface 



We conceive a microscopic model for this setup and describe the superconductors 
by square lattices, i.e., © N^, where the distance between the superconductors is of 
the order of the lattice site distance. The interactions in this system are given by local 
Hamiltonians Hj^, i.e., for a finite lattice with N"^ + N"^ points, 

Hn = Hi^N + Hii^N + Vat. (1) 
Here, Hj jsi and Hjj p^ are the Hamiltonians of the two superconductors, labelled by / 
and //, and Vn models the interaction of the junction between the two superconductors. 
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It is important to realise that the individual superconducting plates are two dimensional 
and that the contact is onedimensional, but nevertheless infinite in the thermodynamic 
limit, i.e., the limit oo. 

The superconductors are modelled by the strong-coupling BCS model, using the 
quasi-spin formalism fOlEI, i-e., the interactions in the superconductors are given by 
the Hamiltonians i?/,iv and Hjj i\f for finite iV e N, 

N-l ^ N~l 

H,,N = Yl ^i^tik^l) -jf^H at{k,l)ar{m,,n,), e, > 0, ^ = /, /J. (2) 

fc,Z=0 k,l,m,n=0 

These Hamiltonians Hj^ act on the Hilbert space (S)j=o^'' (^t^ '^i ^^'^ '^t — 
(jf(j^ — (y^uf are copies of the well-known Pauli matrices, generators of the algebra of 
2x2 complex matrices M2. The operators af{k,l) and aj{k,l) represent respectively 
the creation and annihilation operators for a Cooper pair in superconductor / with 
indices (/c, Z); the are the kinetic energies of the Cooper pairs. 

The junction between the two superconductors / and // is modelled by 

N-l 

^^ = -^ H cri{h,0)ayj{k2,0) + aj{h,0)aij{h,0), 7 G M+. (3) 

fcl,fc2=0 

It describes the rate by which Cooper pairs tunnel through the barrier (Fig. Q). A 
Cooper pair with indices (fci,0) in the first superconductor can tunnel through the 
junction and become a Cooper pair with indices (^2? 0) in the second superconductor. 
The coupling constant 7 governs the rate at which these processes occur. Note that 
only Cooper pairs which are on the contact surface, i.e., pairs which have a second 
index equal to zero, can tunnel through the junction. The BCS model Q we use here 
is permutation invariant with respect to the lattice indices. The underlying geometry 
of the lattices plays no further role. The important point is that only N out of N'^ sites 
of one superconductor, see equation 0, are linked with the other superconductor. The 
high degree of permutation symmetry in the system makes the model exactly solvable 
in the thermodynamic limit, N —>■ 00. 

2.1. Equilibrium properties of the noninteracting superconductors 

Before studying the coupling between the two superconductors 0, we briefly 
describe the equilibrium properties and the phase transition in the noninteracting 
superconductors, i.e., the situation 7 = 0. We treat the first superconductor explicitly, 
the second is analogous. 

Exploiting the permutation invariance of the BCS Hamiltonians (j21) with respect to 
the lattice indices, the properties of the equilibrium states in the thermodynamic limit, 
— i> 00, of this model can elegantly be derived using ^21- The extremal /5/-KMS or 
equilibrium states at inverse temperature /5/ in the thermodynamic limit are given by 
the product states a;^^, where the expectation value oi aW. X = Xq ® Xi ® X2 ® ■ ■ ■ vo- 
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the infinite tensor product algebra of local observables, Bj = [Jbcn^ '^ieB ^2,?., is given 
by: 

^4>iix) = n ^"^pc^iA^i- (4) 

Here, p^^^j is the i-th copy of the 2x2 density matrix p^^ G M2, a solution of the 
selfconsistency or gap equation: 

exp[-^//i^,] , . 

where /i^^ is a one-site effective or selfconsistent Hamiltonian, 

h^, = eja} - Xie'^'aj - A^e-^-^^a/, (6) 
with an order parameter A/, satisfying the selfconsistency equation: 

A/ = |trp<^,a/| = \uj^j{aj) \ ; (7) 
derived from equations (jl} and and with a phase 0/: 

0/ = argtrp<^^cr/ = aiguj^j{aj), (pi e [0, 27r]. (8) 
By explicit calculation, the gap equation ((Zj) can be rewritten in the form 

A,(l--^tanh[/?,fc,]) =0, (9) 



where kj = ^yej + Xj and {—kj,kj} is the spectrum of the effective Hamiltonians h^j 
(jni), and is independent of the phase 0/. It can readily be seen that the equation Q has 
always a solution A/ = 0. It yields the normal phase state. Solutions A/ 7^ exist if the 
following conditions are satisfied: 

ei < 1/2, 

(10) 



These solutions Xj 7^ 0, 0/ G [0, 27i] are understood to describe the superconducting 
phase states. The order parameter A/ 7^ ((Tj) depends on the inverse temperature 
f3i which has to be larger than the critical inverse temperature for the BCS model (3c, 
cf. dim). The phase of the superconductor 0/ ((HI) can be fixed freely. This leads to an 
infinite degeneracy of the states in the superconductive regime and is due to spontaneous 
symmetry breaking fHlElE]- The second superconductor II has analogous normal 
and superconducting equilibrium states. 

2.2. The nonequilihrium steady state (NESS) 

Suppose that initially the two superconductors are both in a superconducting pure 
phase equihbrium state at inverse temperature /5/ and phase 0/ for the first BCS 
superconductor and inverse temperature and phase 0// for BCS superconductor 
//. I.e., the state of this system is the product of two pure phase states ((H), ui^^ ® oJ^jj^^ 
on the total algebra of local observables B = IJ_BcN2eN2 ^ieB 
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If the interactions between the two superconductors are switched on, the system 
should ultimately evolve to a nonequilibrium steady state uj for the dynamics of the 
total interacting system ((H), i.e., 

UJ = lim uj^^ (g) u^jj o at (11) 

t— >oo 

where {at} is the dynamical semigroup of the total interacting infinite system, i.e., 

at{.) = lim e'*^^. e"'*^'^. (12) 

However, it is not clear that the limit (fTT|) exists for our model (^. Computer simulations 
suggest that the system continues to evolve periodically in time. Furthermore, it has to 
be specified in which sense the limiting dynamics at (I12j) can be defined. 

To deal with these problems, Ruelle proposed the following construction for a 
natural nonequilibrium steady state [SHEI: 



1 



T 



w = lim — / dt uqo at, (13) 







where uo is the initial equilibrium state. Ruelle considered the case of a finite system 
coupled to a number of infinite reservoirs, and the conditions on the interactions between 
the small system and the reservoirs where such that the limiting dynamics (jl2j) exists in 
norm on the C*-algebra, i.e., independent of the state of the system. Clearly, the state 
u) (I13p is an invariant state for the dynamics of the total interacting system (|12p and it 
is called a natural nonequilibrium steady state. In recent years, exact results concerning 
entropy production in quantum statistical models were derived using Ruelle's framework 
dg, see e.g., 011111113 • 

In our model (^, the interactions are not local and the existence and properties 
of the dynamics (jl2j) depend on the state of the system. Hence, another approach has 
to be followed to obtain a nonequilibrium steady state. For this model, we propose 
to construct a steady state by selecting a good state which satisfies all the essential 
properties a NESS should have. Which are those properties? 

First, due to the high degree in permutation symmetry in the model (0), it is natural 
to look for steady states in the class of product states, i.e., 

Cd{X)= Yl trp.Xi, X = Xk®Xi®Xm®---eB-XjeM2. (14) 

Furthermore, we expect the state uj to have the same permutation symmetries as the 
total Hamiltonian (^. Analysing this, we see that the total system should be divided 
into four parts (see Figure |21) which are permutation invariant subspaces. The first part 
is the bulk term of the first superconductor (/„), the second piece is the contact or surface 
part of the first superconductor (/;,), the third part is the contact part of the second 
superconductor {Ih) and the fourth part is the bulk part of the second superconductor 
(J/q). The size of the bulk subsystems (/„, //„) is of the order of 0{N{N — 1)) and the 
subsystems on the contact surface (/ft,//^) are of the order of 0{N). Hence, the states 
uj are products of four symmetric product states on these different regions: 

UJ = uJi^® uJi, (g) uJn^ (g) uJn^ . 
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Figure 2. Devision of the system into four permutation invariant subsystems 



hi 



Secondly, requiring that is a steady state is mathematically translated into the 
following expression: 

lim cD([i/7v,X]) = 0, \/XeB. 

AT— >oo 

Due to the product structure of the steady state (fT^ . the limit ^ oo exists; it yields 
an equivalent Hamiltonian of the type Hn = Yli ^ii the condition for the time 
invariance becomes now: 

Cj{[huXi]) = Q, VXiGMs,,, (15) 

where G M2^i is the site dependent effective one-site Hamiltonian for site i G N^©N^. 
It follows directly from the permutation symmetries of (0) and the product structure 
of the state uj (jl4j) that these effective Hamiltonians have the following structure: 

e/a|(i) - u{ajjaj^{i) - u{aYja+{i), i G 4; 

eialii) - u{afJaY^{i) - u{aYjal{i) 

--f(u{aij )aY {{}+u{aYj )ai{i)), i G h; 

where af is the t-th copy of the Pauh matrix a* in the site t G {/a, /&, //a, Ih}- The 
expectation values ijj{(jf) appearing in these expressions are the order parameters for 
this effective Hamiltonian. They are determined in a selfconsistent way, see below, by 
the effective Hamiltonian h (fT^ and the invariant state uj (fT3j) . 

A third requirement is that a NESS uj should have the same properties as the initial 
equilibrium state uj^j ® uj^j^ if we look at the system far enough away from the contact 
surface, i.e., the expectation values for local observables X G with a support in or 
Ha satisfy: 

uj{X) = uj(f,j{X), for X with support in la] ^^^^ 
uj{X) = ujfi,jj{X), for X with support in Ila- 

In this sense, the steady state uj has a memory of the initial equilibrium state uj^j ^uj^jj, 
a property which seems indeed a natural ingredient for a NESS. 
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The NESS can now be constructed explicitly as follows. The requirement that lj is 
an invariant product state (fT3j) for the dynamics induced by h fll6p leads to the following 
choice for the one-site density matrices pi of uj ()14|1 : 

= c+|+)(+|i + c_|-)(-|i, Vi G N^eN^, 

where and |— )j are the eigenvectors of hi, and the constants c_ and c+ are two 
positive constants adding up to one. A natural way to implement condition (fTTj) is to 
to define the constants c_ and c_|_ by means of the expectation values of the one-site 
density matrices of the initial equilibrium state tu^^- ® on the eigenstates of hi, i.e., 
we chose 

c+ = {+\po,i\+)i, and, c_ = {-\po,i\-)i- 

In other words, for every site i in both superconductors, we project the local one-site 
density matrices po,t of the initial equilibrium state uj^j- ®uj^jj onto the subalgebra in M2 
which is invariant under the dynamics induced by h Hence, the one-site density 

matrices pi of the NESS are given by 

A = I+)(+I(+,Po+)^ + |-)(-|(-,Po-)m VzGN^eN^. (18) 

Furthermore, in order to meet condition (fTTj) and to avoid trivial nonsuperconducting 
solutions in (fTH|) . we explicitly set the bulk order parameters and phases equal to the 
equilibrium order parameters and phases, i.e., 

(:;(a+) = A,e'<^^=^<^,(a+), and, ^{aj.J = Xje't'" = uj^,,{aj,). (19) 

It is easy to check that this choice ()18 p - (fT^ yields a state uj fulfilling indeed all of 
the above requirements. By construction, this state is again a product state (fT4|l and it 
is an invariant state under the limiting dynamics (fT3|) . In sites which do not belong to 
the contact surface, i.e., i ^ (^i, 0) G h and i 7^ {k2, 0) G the local density matrices 
of the initial equilibrium state and the constructed NESS coincide (fT7j) . This can be 
seen as follows. In the bulk regions, i.e., for sites i G laUlIa, the effective Hamiltonians 
hi (fTHjl coincide with the equilibrium effective Hamiltonians h^pj or /i^^^ (jH)). Hence in 
la U Ha, the spectral decompositions of p^j or p^^^ coincide with those given in (fTHj) . 

Now, we have to compute the order parameters and phases in the effective 
Hamiltonians Compared to the equilibrium situation, see section where two 

order parameters ((Tj) and two phases ^ completely determine the equilibrium states, 
we have now four different order parameters and four phases appearing, cf. (jlf)|l . i.e., 
apart from the parameters in Equation (fT^ . we also have to compute 

Xi = |a'(o-/J|, and, 0/ = arga;(o-/J; 
Xjj = \uj{(rfjj\, and, (pu = aTgCj{afjJ, 

Using ()18j) and (fTIH) . we find an explicit expression for the selfconsistency equations in 
the steady state u) which determine these order parameters and phases: 

Xje'^^ =uiai{h,0)) 
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^ef + A/|A/e'<^^ + 7A7/e*<^" | cos 



(A,e^<^^ + ^\ne'^- r ^ "^^^ , ^ ^^^^^^ , ' ^^Z ~ (20) 



= (A..e'^- + ^j^^,i^.)^?. + A..|A..e^^^^+7A.e'^-|cos(0..-0..) 

ef/ + \Xiie''t"! +7A/e'<^^|2 

Together with the selfconsistency equations for A/ and Xu the above equations form 
a set of six coupled transcendental equations whose solutions completely determine the 
state uj (HH). 

The steady state u) (fTH|) divides the system into four parts. In both superconductors 
away from the contact surface, see Fig. |2l the presence of the other superconductor 
is not felt. The system is unperturbed and behaves as in the initial equilibrium 
state. Nevertheless, the NESS (fTHjl is not normal with respect to the equilibrium state 
® u^jj , because the one-site density matrices pi on an infinite amount of sites on 
the contact surface (i.e., the sites i E Ih) differ from the one-site density matrices 
of the equilibrium state (0). On the contact surface, the system is effectively disturbed 
and there appear two intermediate layers with different properties compared with the 
initial equilibrium state. 

Finally, let us compare this construction of a NESS with the construction proposed 
by Ruelle (fT!?|) . The essential relation between our construction (fTH|l and Ruelle's 
construction (fT^ is that in the latter one, it can be expected that the parts of the 
algebra which are not invariant under the dynamics ()12|1 are averaged out and that the 
NESS is effectively a projection of the initial state on the subalgebra which is invariant 
under the dynamics, i.e., it leads to the same state as in (fT8| . 

This finishes the construction of our framework. Let us summarise the situation 
at hand. We start with the microscopic model (Q) of two different superconductors 
interacting via a particular Josephson junction and we constructed a particular 
reasonable nonequilibrium steady state Co with site local density matrices given by 
Equations (|TR|l - (jTni . 



3. Transport phenomena 

In the framework described above, we can investigate the nonequilibrium properties and 
the existence of currents which emerge by bringing two superconductors in contact with 
each other. In particular, we compute the expression for the Josephson particle current. 

Suppose Qiv is an extensive observable, e.g., the total number of particles, the 
energy,. . . The transport related to this quantity i{Q) in the steady state (fT8|l is defined 
as, 

j(g)= lim il(:;([i^^,Q;v]), (22) 

where is the Hamiltonian of the total interacting system (P). In this definition, the 
current is scaled with the size of the contact surface, i.e., we divide by A^. This is the 
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right scaling to obtain density quantities. Therefore, j{Q) is the current density. Taking 
into account the product character of the steady state (fTHjl and the fact that this state 
is invariant only under the limiting dynamics we compute 

^ N-l N-l 



fci,/i=0 



k2,h=0 



N{N-1) 



N-l 



N-l 



J2 (ri{h,h)\je-"'" + aY{h,h)\je"^' 



ki,h=0 



'N S (Tt{ki,li)~Xie-"^' +aY{h,h)~Xje"^ 



fci,/i=0 
N-l 



fci=0 



N{N-1) 



N-l 



Ar2 



N-l 



k2,l2=0 



'N S ^tiik2,l2)Xiie-"^" +cTij{k2,l2)Xiie''^ 

k2,h=0 
N-l 

fc2=0 

This is simplified using that uj is invariant under the limiting dynamics, i.e. 
^{[T.iKQN]) = (US). We obtain 

N-l N-l 



ki,li=0 k2,l2=0 



N 



(23) 



N 
where 

Xj{h,li) = (ri{h,h)CXie'^' - A/e'-^O + (rj {h,h)CXie-'^' - A/e'^^O, 
and 

Xn{k2, 12) = crUk2, /2)(A//e'^" - A,,e'<^") + ajjik2, l2){Xue~'^" - Xue-"^")- 

The order parameters and phases. A/, A/, . . ., appearing in these expressions are defined 
in Equations and (j^U- Already from this expression (j^ . it is clear that the 

general current density j{Q) for arbitrary Q (see Eq. (12 is of order 0(1). In the 
following sections, we compute the transport of two important quantities. 



3.1. Particle current 



Of course, the most interesting transport phenomenon in this model is the Cooper pair 
current through the junction. Is their really a current and what are its properties? The 
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J{C) 




Figure 3. Dependence of the Josephson current j{C) on the phase difference A<j> 



relative particle current is defined as the transport of the relative number operator Cn 

N-l N-l 

Cat = ^ aj{ki,li)aY{ki,li) - ^ cr/f(fc2, ^2)0-7/(^2, ^2) 

ki,h=0 fc2,«2=0 

Calculating j(C) using equation (j^ . we find 

j{C) = 2{Xi~Xi sm{<Pi - 4>i) - XiAii sm{<j)ji - 4>ii)) (24) 

The order of magnitude of the total flux of Cooper pairs is proportional to the size of the 
contact surface. It depends explicitly on the sines of the phase differences between the 
bulk and the contact surface, i.e., (pi — (pi in both superconductors, i = I, II, and on the 
order parameters, see (fT^ - (PT|) . However, the only free parameters are the bulk phases 
(pi and (pi I, the other parameters are determined by the selfconsistency equations pOjl - 
()21|). Hence, the interesting point here is the dependence of the Cooper pair current 
on the phase difference between the two (bulk) superconductors A(p = (pu — (pi. Let 
us illustrate this with an example. In Figure El we plotted the Josephson current j{C) 
as a function of the phase difference A0 = (pu — (pi. For this figure, we took 
/?/ = I3ii = 4, ei = en = 1/4 and 7 = 1/2. If the phase difference is zero A(p = 0, 
there is no Cooper pair current in the system, and if the sign of the phase difference 
is inverted, also the direction of the superconducting current is reversed. Hence, this 
supercurrent j(C) has the basic properties of the superconducting tunnel current 
predicted by Josephson ;1]. In the textbook derivations of the Josephson current, the 
current is proportional to the sine of the phase difference A(p. In this model, however, 
the relation with the phase difference is more complex ()24|) . but as can be seen in the 
example (Figure Ej), it has a similar sinusoidal dependence on the phase difference of 
the order parameters in both superconductors. 
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3. 2. Heat fluxes and entropy production 

Here, we want to make a link with the current research about entropy production in 
quantum statistical models [H|l7|IT!)j . We investigate the entropy production in our model 
and compare different notions. 

In quantum systems when the time evolution is given by a unitary conservative 
dynamics, it is known that the von Neumann entropy does not change. This unitary 
time evolution however, induces correlations between different parts of our large system. 
It is natural to think of an observer who can only observe part of the large system. When 
after some time the correlations of this part with other parts in the large system are 
forgotten, the entropy of this subsystem has changed, or entropy has been produced. 
In this way, entropy production can be understood as the change of entropy when 
partitioning the total system in different parts [T?)] . 

Applying this idea to our model and considering the entropy production induced by 
the partitioning of the total system into the two superconductors, we find the following 
expression for the entropy production 16J: 

d = lim -i^u{[HN,logij\j + \oguj\jj]), (25) 

N^oo iV 

where u\j and are the density matrices of the state u restricted to resp. the first 
and the second superconductor. 

However, due to the product character of the steady state and using cyclic 
permutations under the trace, it is straightforward to check that the expectation value 
of this commutator vanishes. Hence, there is no entropy production in this model, i.e., 
ei = 0. 

The definition can be compared with the conventional notion of entropy 
production, e.g., see [SIIZ|, where the entropy production is defined as the sum of the 
heat fluxes in the system, i.e., we calculate the current ()22|1 associated with the following 
observable 

where Pi and Pu are the initial inverse temperatures of the equilibrium states of the two 
superconductors and Hj^n and Hjj n are the BCS Hamiltonians (j2I). They represent 
the internal energy in the superconductors. This leads to the following current ()22|) : 

62 = j{Q) = lim ^3j{u{a'j{h, 0)) + 2ej)Xi~Xi sm{(j)j - ^j) 

+^/5,,(cu(a|,(A;2, 0)) + 2en)Xnhi sm{^n - hi) (26) 
= 0. 

Hence, the first definition ()25|) of entropy production and the conventional definition 
fl22|) in terms of heat fluxes both lead to zero entropy production. However, for the last 
case, we see that the heat fluctuations are of order 0{N~^) around zero. Also, there 
is no definite sign for the entropy production fluctuations in this state as the sign 
depends on the phase differences, i.e., just as in the case of the Josephson supercurrent 
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(|24p ■ but these heat currents are of a lower order of magnitude. This way, we may 
conclude that there are no heat fluxes in this system, there are only small fluctuations 
around zero. 

Finally, remark that in our setup, we took two dimensional superconductors and a 
onedimensional contact surface. Of course, it might be interesting to look at other 
dimensions for the superconductors and the contact surface. We expect that our 
methods and results remain valid in those situations. The onedimensional case, however, 
is basically different. If the contact surface is finite in the thermodynamic limit, the 
nonequilibrium steady state does not factorise in the contact region. 
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